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Abstract. In this paper we consider the general setting for constructing Action Prin- 
ciples for three-dimensional first order autonomous equations. We present the results 
for some integrable and non-integrable cases of the Lotka-Volterra equation, and we 
show Lagrangian descriptions which are valid for systems satisfying Shil'nikov criteria 
on the existence of strange attractors, though chaotic behavior or homiclinic orbits 
have not been verified up to now. The Euler-Lagrange equations we get for these 
systems usually present "time reparameterization" symmetry, though other kinds of 
invariance may be found according to the kernel of the associated symplectic 2-form. 
The formulation of a Hamiltonian structure (Poisson brackets and Hamiltonians) for 
these systems from the Lagrangian viewpoint leads to a method of finding new con- 
stants of the motion starting from known ones, which is applied to some systems found 
in the literature known to possess a constant of the motion, to find the other and thus 
showing their integrability. In particular, we show that the so-called ABC system is 
completely integrable if it possesses one constant of the motion. 
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1. Introduction 

Three-dimensional autonomous first order nonlinear systems show a very rich behavior, 
from complete integrability [I], |3], [| to chaos and strange attractors^], H, [7L 
according to the values of the parameters that appear in the equations of motion. 
As they usually represent physical and biological systems of interest (Lorenz equation 
for hydrodynamic flow, Rossler equation for chemical reactive systems, Lotka-Volterra 
equation for laser physics, plasma physics, [|T0|| biology, ]ll], Q economy, Hl3| etc.), it is 



somewhat surprising the lack of Action Principles for dealing with these equations. We 
intend to fill this gap using the general setting for Lagrangian structures, which 
when applied to three-dimensional systems reduces to a simpler problem, namely: for a 
given three-dimensional vector field V 8 ", to find the determinant of a metric (or volume 
element) such that the covariant divergence V a ^ a = 0, which is a way of writing 
Liouville's theorem and also a condition which is known to lead to the construction 



of constants of the motion starting from symmetries only|T5|. The Euler-Lagrange 



equations that come out from these principles are usually equivalent to =j- = V a , where 
s is an arbitrary parameter (time reparameterization symmetry). 

In section |2| a preview of known and new results on the Lagrangian approach for 
three-dimensional systems is presented. Section |3| deals with applications of the theory: 
a method to construct new constants of the motion from known ones, construction of 
Hamiltonian structures, and several kinds of Action Principles with invariance related to 
symmetries. In section |] we show examples for the Lotka-Volterra equation, including 
integrable cases as well as non-integrable ones: we show that recently found "quasi- 
integrable" cases|| (i.e., systems possessing one constant of the motion) of Lotka- 
Volterra are indeed integrable; for integrable cases, the construction of Lagrangians 
leads to new related equations, involving symmetries (of the evolution vector V) which 
generate the kernel of the symplectic 2-form. The Euler-Lagrange equations that arise 
possess several types of invariance, the simplest of which is time reparameterization, 
which arises when the kernel of the symplectic 2-form is generated by the evolution 
vector V; finally, two examples (related to known chaotic systems in the literature^, |5|) 
show that Lagrangian descriptions for systems near chaotic behavior do exist and they 
satisfy Shil'nikov conditions |7|] on the existence of strange attractors around the fixed 
points, though it is still uncertain whether a homoclinic curve exists or not for such 
systems. 

2. Preview 

Consider the autonomous equations of motion: 

x a (t) = V & [x h {t)} , a, b = 1, 2, 3 , (1) 

where V defines a flow, velocity or evolution vector field in the corresponding space, 
being a function (which is local in time) of the coordinates x a (t). As an example, the 
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Lorenz equation is defined by the vector 

V[x, y, z] = {V x , V y ,V z } = {—a x + ay, r x — y — x z, xy — bz} , 
where (a, r, b) are three independent parameters. 

2.1. Lagrangian structures for three-dimensional systems 

In arbitrary dimensions, the problem of finding Lagrangian formulations or Action 
Principles for the given equations of motion (|IJ) leads to covariant equations [14] in terms 



of Lie derivatives of tensor structures along the flux vector. In three dimensions, the 
equations are very restricted: as it is an odd-dimensional system, the Euler-Lagrange 
equations always have an invariance. For example, time reparameterization invariance 
appears when the kernel of the symplectic 2-form is generated by the evolution vector 
field V. 

From now on, unless explicitly stated, we will be dealing with objects which are 
time— independent, in the sense that they only depend on time implicitly, through 
their dependence on the coordinates x a (t) . 

The equations of motion (|I]) are related to a Variational Principle with Action 

S[x a {t), t] = ^ dt (L a (x a - V a ) + K) (2) 



(we assume from here on a summation over repeated indices), where the 1-form L[x h ] 
and the zero- form ^[V 3 ] satisfy the following equation: 



L ath V h + L h V\ a = K, a , 



with K, a ee g . 

We rewrite the above equation in terms of invariant structures: 

CL = dK , 

v 

where C is the Lie derivative along the vector V, and d is the exterior differential 

v 

(see [|T6|l for a definition of these operators). 

We call the pair (L; K), 1-form L and 0-form K, a standard Lagrangian pair 
for V if K ^ 0. In the special case K = we call L a non-standard Lagrangian for 
V: this case allows for the construction of Poisson brackets and constants of the motion 
(see section ||). 

The Euler— Lagrange equations which come from the Action (0) are: 

s ab (x b -y b ) = 0, (3) 

where £ = dL is the symplectic 2— form or Lagrange bracket whose components are: 

Notice that these Euler-Lagrange equations do not imply the original equations of 
motion (]!]); instead they imply deformed or mixed equations involving the kernel of 
the symplectic 2-form. In three dimensions, the symplectic 2-form always has a kernel, 
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and the simplest thing to do is to ask the kernel to be proportional to the evolution 
vector V itself (for other types of kernel see section p.3| ): 

Sab = p[x d ]e ahc V c , 

where e a t, c is the Levi-Civita 3-form, a totally antisymmetric tensorial density, with 
e i23 — 1 j and p[x a ] is the volume element, a scalar density which solves |H)| 



(Pn,a = 0, (4) 

which is sufficient and necessary in order that S be a closed symplectic 2-form for the 
evolution vector V. In terms of a metric g a b, the volume element is p = \J det g a b and 
the above equation is equivalent to V a ^ a = 0, where the covariant derivative is taken 
with respect to the metric g. Geometrically, the existence of such metric leads to the 
construction of an integral invariant [|17| (Liouville theorem) of the form 

p[x a ] d 3 x , 

M 

where M is a region of the configuration space which is assumed to evolve along the flux 
V. The volume element is not unique: in fact, given two volume elements p,p', their 
quotient C = p/p' is a constant of the motion. 

As a consequence of the fact that the kernel of the symplectic 2-form is generated 
by the evolution vector, our Action Principles for three-dimensional systems will be 
purely geometrical: they will give only the curve of the motion, leaving the way local 
clocks run undetermined. The Euler-Lagrange equations (H) are thus S a bi b = , which 
are equivalent to ^f- = V a , where s is an arbitrary parameter. 

The standard Lagrangian pair (L; K) is constructed (locally) by a line integral: 

K[x b ] = L a [x h ] V a [x h ] , 

L4x b ] = fix ah [ x *( s )}^sds + R a [x*}, U 

where R is an arbitrary 0-form and the path in the line integral is parameterized by 

x a {s) = x a + s {x a - x a ) , s G [0, 1] , 

and is such that the integrand is equal to zero at s = and well behaved along the path. 
The Action for the system shows clearly time reparameterization invariance: 

S[x a (t)}= L a [x h (t)]x*(t)dt. 

It is important to mention here the closeness of the equation for the volume element 
with that for the Darboux polynomial. In the usual terminology [0, two polynomials 
in the coordinates f(x a ) (the Darboux polynomial) and Q(x a ) (the cofactor) are to be 
found in such a way that 

f,,V a = Qf, 

therefore we may get a solution for the volume element, equation @ from a Darboux 
polynomial when the cofactor is a multiple of V a ^ a . There are usually other solutions 
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for the volume element such that the resulting expression for / is non-polynomial, as it 
happens with the so-called exponential factors in Ref. 

Anyway, this closeness may be used directly to construct Action Principles for 
every equation for which a Darboux polynomial or "Darboux function" is known, as it 
happens in the Lorenz equation, the Rossler equation, the Lotka-Volterra equation and 
many others. 

3. Applications 

In the following applications, we assume that the evolution vector V possesses a non- 
standard Lagrangian L (this fixes R in equation (H)). 

3.1. Method to construct a constant of the motion starting from a known one 

We may use a known constant of the motion to construct the other explicitly. We will 
depict the method, omitting the details for simplicity. 

Assume that the evolution vector V possesses: (a) A non-standard Lagrangian L, 
along with its associated volume element p ; (b) a time-independent constant of the 
motion if[x a ] . Consider the vector defined by the components 

v a = - e abc if h L c . 
P 

It follows that r) is a symmetry for V, and that rj is proportional to V. Therefore 
we obtain 

rf = I V a , 

where J is a constant of the motion for V, which may or may not be independent 
of H. We will show that the system is completely integrable (i.e., it possesses two 
independent, time-independent constants of the motion) in any case: 

(i) if i = 0, then DH a = L a where D is an if -independent constant of the motion, 
because, otherwise, dL would vanish, 

(ii) if I 7^ is independent of H, it is direct and 

(iii) if I ^ is a function of H, I = 1(H), define the 0-form P(H) s.t. % = j^) ■ 
Define next the 1-form U = e~ p L . It follows that dll = 0, therefore the 0-form C 
defined (locally) by 

dC = e- p L 

is an if -independent constant of the motion for V. 

3.2. Hamiltonian theories 

Using the constants of the motion obtained in the last subsection, it is easy to construct 
all possible Hamiltonian theories for the evolution vector V. For example, for cases (ii) 
and (iii) above we may write 
V a = J3 b H b ^ 
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where J ab = e abc L c is the Poisson Bracket. The Casimir, though always 

computable, is easy to guess only in case (iii): it is the above 0-form C, and it solves 
</ ab C,b = 0. The Hamiltonian is H modulo a function of C; its final selection relies 
on stability conditions (see Ref. p8|). 



3.3. More Lagrangian Theories 



If we have two independent constants of the motion for system ([!]), then another time- 
dependent constant of the motion may be constructed (in principle) by integration of 
any component of the evolution vector. If H l [x,y,z], H 2 [x,y,z] are two constants, 



assume we may solve for y and z, getting y = y[x, H 1 , H 2 



z[x,H l ,H 2 ]. Thus 



the first components of the equations of motion, x 
x = Fix,!! 1 , H 2 } and thus 

H 3 [x,y,z,t] = H 3 [x,y,z}-t 



V l [x,y,z\, are transformed into 



is a time-dependent constant of the motion, where H 3 [x, y, z] = y J x (F[x, H 1 , H 2 ]) dx 
Now that we have three constants of the motion, we assert that the following is the 



generic form for any non-standard Lagrangian 1-form: |L4 

L = C 1 dC 2 , 

where C j = C j [H 1 , H 2 ,H 3 -t], j = 1,2. 

In the case that any C J depends on if 3 , the associated kernel of the symplectic 
2-form is not anymore proportional to V: indeed it is proportional to a (probably 
time-dependent) symmetry rj for V, which may be calculated explicitly. We get, as a 
result, new Euler-Lagrange equations which mix the vector V with rj. Examples with 



time-independent symmetries r\ will be shown in section 4.2 



Hi=Hi[x,y,z\ 



3.4- Symmetries 

Finally, assuming that the 0-forms H l , H 2 , H 3 from the last subsection are known, then 
an Abelian algebra of vectors A = {r]i, r] 2 , r] 3 = V} is constructed by taking the dual 
vectors of the 1-forms dH l , dH 2 , dH 3 , i.e., by solving for the vectors: 

CH k = 5 k 3 , j, k = 1,2, 3. 

These commuting vectors are symmetries of each other by definition, and they are 
used to generate the kernel of the most general symplectic 2-form for the evolution 
vector V (see section p~2| ). 

4. Examples: The Lotka— Volterra equation 

This equation, restricted to three dimensions, reads in its general form, 

x = V x [x,y,z] = x (ai + b n x + b 12 y + b 13 z) 
y = V y [x,y,z) = y{a 2 + b 21 x + b 22 y + b 2Z z) 
z = V z [x,y,z] = z (a 3 + b 31 x + b 32 y + b 33 z) , 
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where constant parameters. We look for a volume element of the form 

p[x,y,z]=x u y v z w , (6) 

where u,v,w are constants that depend on the above parameters. 

For some values of the parameters, non-standard Lagrangian 1-forms may be found: 
their existence allow to define Poisson brackets, and show explicit integrability. In 
other cases, only standard Lagrangian pairs are obtained, though the question if these 
Lagrangians could be made non-standard by an addition of a closed 1-form is open. The 
examples suggest that, for chaotic or near chaotic systems, the answer to this question 
is negative. 

4-1. Case ba = , i = 1, 2, 3 (no Verhulst terms) 

The evolution vector is, after rescaling, 

V[x,y, z] = {x (A + C y + z), y (fx + A z + x), z (v + B x + y)} , (7) 

where A, B, C, A, ft and v are constant parameters. Constants of the motion for this 
system exist for a subset of the parameter space. 0, ^ |], U For arbitrary values of the 
parameters, however, this system is not integrable, but as we will see, a Lagrangian 
description always exists. Consider the volume element 

p[ x a ] = (xyz)- 1 . 

It is easy to see that (p^ a ), a = , and the Lagrangian 1-form is obtained directly: 

L a [x h ] = (xyz)- 1 e ahc (x h V c [x d ] + V h W c [x A ]) : (8) 

where if = M a b x b , M a b = diag(A, fx, v) , and W 3, = x a lnx a . 

Notice that the above volume element is singular at the planes x = 0, y = 0, z = 0, 
reflecting the fact that these are invariant planes. 

4.2. ABC system 

When A = fx = v — , the system defined by the flow (0) is called "ABC system" . In 
this case, the Lagrangian (|[) is non-standard, and reduces to 

L a [x h ] = (xyz)- l e lihc x h V c [x d ]. 

The fact that this Lagrangian is non-standard leads us to conclude: if (for given 
values of the parameters) there is a constant of the motion, then another may be found 
using the theory of section |3|. This may be applied, for example, to all the ABC systems 
studied in Ref. |3J where the author finds one polynomial constant of the motion. In the 
next examples, the method will be applied to some new cases found in Ref.||. 

As a final aside, the Lagrangian may be used to find Poisson Brackets and 
Hamiltonian theories, and to find new, related equations of motion from the construction 
of a wide class of Lagrangian theories, as it will be done in the last example of this 
subsection. 
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(i) Case A = — 1, B = 1/2, (7 = 0. 

This case possesses a recently found constant of the motion: H 1 = 
xy~ 1 z 2 exp (— 2(y + zY^xy)^ 1 ). We apply our Lagrangian to construct another 
constant of the motion. In the present case, it turns out that U = (H 1 )^ 1 ^ 2 L 
is a closed 1-form, which therefore must be the exterior derivative of some constant 
of the motion. We obtain, after integration, 

r(y+z)/y/xy 

H 2 [x, y, z] = (H l Y l/2 x-2 exp (q 2 ) dq , 

Jo 

and thus this system is integrable. 

(ii) Case A = -1/(C + 1), B = 2. 

This is a new case also found in Ref.||, with one constant of the motion: 
H l [x, y, z] = x 2 \y\ 2 ( c+lS) \ z\ ~ 2C \ 2A 2 x z — (y — A z) 2 \ c ~ l . According to our theory, in 
this case we obtain that the 1-form U = (H 1 )" 2 <. 1+c '> L is closed. We integrate and 
get the other constant of the motion: 

H 2 [x, y, z] = (y - A z) (l - n + x y tl AC ( 1 - fi) ^ A ^^c) ®n (-AC, \- A^j 

where = and Bn(a, b) = ^ q a ~ l (1 — q)^ 1 dq is the incomplete Beta 

function. 

(hi) Case ABC -1 = 0, B(A + 1) + 1 = 0. 

Here a quadratic constant of the motion[|T], |3| is known: H 1 [x,y,z] = A 2 (Bx — 



z 



,2 



2 A(B x + z)y + y 2 . We obtain the other constant: 

1+2 C 

H 2 [x,y,z] = \x\~ 1 \yz\~ 1 ~ c x ±\ r W [ + A(B x - z) - y 



x 



(y + A z)x - (y - A z)C(±VW T + y-Az) 



The above three systems are, therefore, integrable and the construction of many 
Hamiltonian as well as Lagrangian theories is possible. In the next example we will 
take a known integrable case to show the construction of Lagrangian theories which 
give new equations of motion, in terms of the symmetries of the evolution vector 
field. 

(iv) Case ABC+ 1 = 0, (7=1. 

Here, the integrability is guaranteed by the first condition ||. We use the restriction 
(7 = 1 here for simplicity only. Known constants are H 1 [x, y, z] = — x + y — Az and 
H 2 [x, y,z] = xy B z~ l . Now we take the ^-component of the equation of motion 
(HI), to obtain the third (time-dependent) constant of the motion: 

H 3 [x,y,z,t] = H 3 [x,y,z}-t. 

where 

In ( v ) 

H 3 [x,y,z} = 
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The vectors 771,772,^3 dual to the 1-forms dH 1 ,dH 2 ,dH 3 are commuting vectors, 
and they are found to be, after rearranging: 

H l rfr[x,y,z] = {x,y, z] + H 3 r) 3 [x,y, z] - B H 2 r} 2 [x,y, z] , 
r) 2 [x,y,z] = {x + Az)- 1 y- B z 2 {A,0,-l}, 
V3[x,y,z) = V[x,y,z] = {x(y + z), y(x + Az), z(B x + y)} . 

It is clear from the definitions that any of the three equations of motion (labelled 
by j) which the above vectors define, 

x a = 77/[x b ], J = 1,2,3, 

are completely integrable, with three constants of motion given by Cj = H k — 
t5 k , k = 1,2,3. 

Now we turn to the Lagrangian descriptions for the evolution vector 773 = V (the 
same could be done for the other vectors). Among all the possible examples, 
three special Action Principles are obtained with the Lagrangian 1-forms given 
by L 1 = H 2 dH 3 , L 2 = -H^H 3 and L 3 = H l dH 2 . They are all non-standard 
Lagrangians for V, and the kernel of the symplectic matrix dV is easily found to 
be proportional to the vector r]j. Therefore, the Euler-Lagrange equations coming 
from the Lagrangian U are 

where a is an arbitrary 0-form. The case j = 3 gives the usual time 
reparameterization invariance, but the cases j = 1,2 are examples of other kind of 
invariance; their Action Principles are easily computed from equation (0). 

4-3. Lagrangians and chaos 

In the following examples we construct Lagrangians for Lotka-Volterra systems which 
are close to known chaotic systems. The general results here are: first, that a volume 
element of the form (|6|), implies \^ a a = at the finite fixed point (or singular point) 
of V which is not contained in any coordinate plane. This is somewhat surprising, 
because the volume element is singular at the coordinate planes, but its existence implies 
conditions on the vector field at a point which is far from the planes. 

Second, for a large subset of the parameter space, these systems allow for 
Lagrangian descriptions along with Shil'nikov conditions on the existence of a strange 
attractor of spiral type0 in the vicinity of the relevant fixed point, though no homoclinic 
curve has been found yet. 

Third, the volume element becomes ill-defined (i.e., the powers u, v, w in equation 
([]) become infinite) when the parameter values are such that some fixed points of the 
vector field degenerate into a "fixed line" (i.e., a line in the configuration space for which 
the evolution is frozen) and the system gets a constant of the motion, which may be 
calculated explicitly. 

(i) A replicator-like equation. 
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This example may be understood in the context of catalyst replication and 
mutation: 

x = V*[x,y,z] = x(\(l-x)+\(l-y) + ±(l-z)) 

y = Vy[x,y,z] = y (-1(1 - x )- ±(1 -y) + ±(l- z )) (9) 

z = V z [x,y,z] = z(\x + ti(l-x) + ±(l-y) + ±(l-z)), 

where A is a real parameter, and /i = — ^ — A . The volume element takes the form 



1+6 A 1+6 A 1-2 A 

p = x 6A y 3X z 2A 



It is shown that (p^ a ), a = . The Lagrangian 1-form is found to be 

{L x , L y , L z } = p{ (1 + 5 x — 3y — 3 z) y z + (—30 + 60 x — 9 y) y z A, 

0, 3xy (-11 + 5x + 5y + z) +A5xy 2 \}. (10) 

We note that the case A = and /i arbitrary is found in the literature, displaying a 
one-parameter family of strange attractors||. A bifurcation diagram in terms of \i 
may be found in Ref. 0. We mention it because our case intersects with the latter at 
the point A = , /i = — \ , where the volume element associated to the Lagrangian 
explodes. This point in the parameter space also represents a degeneration of some 
fixed points of the vector field, into the fixed line {1 + 1 — ~, 1 + s}, s 6 M, and 
the following turns out to be a constant of the motion for the system @ : 

C[x,y,z] = \x\ \y\ 2 \z\~ 3 . 

A one predator-two prey system. 

We consider one of the most important models of predation[||] using the Lotka- 
Volterra equation, namely the equation 

x = V x [x,y,z] = x(r - j^x - j^y -bz) 

y = V y [x,y,z] = y(r-j^ax-^y-(b-e)z) 

z = V z [x,y,z] = z(cbx + c(b — e)y — d) . 

The values of the parameters allow to describe competitive superiority of prey x 
over prey y (case a > 1), and predator-avoidance advantage of prey y over prey x 
(case < e < b). Now, if we keep the parameters arbitrary, we get the condition 
for the existence of the volume element: 

o {a — 1) 

and the volume element is 



p i, y, Z] = [X b yb-e Z cK(b-e) 

xyz 



r{a — 1) 



a — a 



where a Q = jr f • 

u b—t b 



In Ref. 0, the author finds a set of values of the parameters for which the system 
develops spiral chaos. [0, || The values are a = 1.5, b = 0.01, c = 0.5, d — 1, r — 
1, K = 1000, and e = 0.009. 

It is easy to check that there is no volume element of the form @ for the above 



set of values, because condition (|ll|) is not met. We get that for the existence of 
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the volume element the only change is d — 8.1, which is much larger that the older 
value, d — 1. According to the model, this means that the predator (z variable) 
has a larger mortality rate when there is a volume element of the form (|^) than in 
the chaotic case. 

On the other hand, a study of this system under condition (|11]) shows that, at 
the fixed point (x ,y ,z ) 6 R_ x E^_, Shil'nikov criteria on the existence of a 
strange attractor of spiral type@] are met, for a large subset of the parameter space. 
However, no homoclinic curve or chaotic behavior have been found numerically yet. 
The volume element becomes ill-defined when a = a Q : again, some fixed points 
of V degenerate into a fixed line and the system gets the following constant of the 
motion: 

C[x,y,z] = \x\-^ {b - e) \y\^ b \z\. 

5. Conclusions 

Lagrangian descriptions for three-dimensional systems are directly related to the 
existence of determinants of metrics such that the covariant divergence of the evolution 
vector is zero, which in turn implies that there is an invariant volume element for the 
system. As a feature of odd-dimensional systems, the Euler-Lagrange equations we get 
usually possess time reparameterization invariance. A Variational Principle for three- 
dimensional evolution equations may thus be useful for the study of (quasi-)periodic 
orbits and long-time properties, or simply to test numerical results. 

In the case of the Lotka-Volterra equation, the very existence of the volume element 
implies some condition on the vector field at the fixed point which is relevant for 
the chaotic attractor. This condition is compatible with Shil'nikov criteria for the 
existence of spiral chaos, but we have not seen chaotic behavior numerically. Anyway, 
the condition does not imply at all the integrability of the system. On the contrary, 
the examples show that the volume element becomes ill-defined when constants of the 
motion appear through degeneracy of fixed points of the flow into fixed lines. 

Recently found quasi-integrable systems (e.g., ABC systems with one constant 
of the motion) are shown to be integrable using a newly devised method to find a 
constant of the motion starting from a known one. Finally, the Lagrangian viewpoint 
for integrable systems leads to Euler-Lagrange equations with several kinds of invariance 
(including time reparameterization), according to the kernel of the associated symplectic 
2-form. 
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